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Abstract. We describe two aspects of the physics of hybrid stars that have a sharp interface between a core of quark 
matter and a mantle of nuclear matter. Firstly, we analyze the mass-radius relation. We describe a generic “Constant- 
Speed-of-Sound” (CSS) parameterization of the quark matter equation of state (EoS), in which the speed of sound is 
independent of density. In terms of the three parameters of the CSS EoS we obtain the phase diagram of possible forms 
of the hybrid star mass-radius relation, and we show how observational constraints on the maximum mass and typical 
radius of neutron stars can be expressed as constraints on the CSS parameters. Secondly, we propose a mechanism for 
the damping of density oscillations, including r-modes, in hybrid stars with a sharp interface. The dissipation arises 
from the periodic conversion between quark matter and nuclear matter induced by the pressure oscillations in the star. 

We find the damping grows nonlinearly with the amplitude of the oscillation and is powerful enough to saturate an 
r-mode at very low saturation amplitude, of order 10“ 10 , which is compatible with currently-available observations of 
neutron star spin frequencies and temperatures. 

PACS. 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions - 26.60.-c Nuclear 
matter aspects of neutron stars - 97.60.Jd Neutron stars 

1 Introduction 

One of the major goals of nuclear physics is to establish the na¬ 
ture and properties of matter at nuclear density and beyond. Ex¬ 
perimental constraints on the properties of ultra-high-density 
matter come from neutron stars, where gravity compresses mat¬ 
ter to the highest densities found in nature. One hypothesis is 
that there is a phase transition from nuclear matter to quark 
matter, and that some stars, called “hybrid stars”, contain cores 
of quark matter. 

We will assume that in hybrid stars the core of quark mat¬ 
ter and the mantle of nuclear matter are separated by a sharp 
interface (Maxwell construction) as opposed to a mixed phase 
(Gibbs construction). This assumption is valid if the nuclear- 
quark phase transition is first order, and the surface tension of 
the interface is high enough. This is a possible scenario, given 
the uncertainties in the value of the surface tension [1-3]. For 
analysis of generic equations of state that continuously inter¬ 
polate between the phases to model mixing or percolation, see 
Refs. [4,5]. 

In this paper we discuss two aspects of the physics of such 
hybrid stars, each of which is relevant to a different set of ob¬ 
servable properties of the star. 

Firstly, in Secs. 2 to 4, we analyze the mass-radius relation 
of hybrid stars, using a a generic parameterization of the quark 
matter equation of state (EoS). We obtain the phase diagram of 
possible forms of the hybrid star mass-radius relation, and we 
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show how observational constraints on the maximum mass and 
typical radius of hybrid stars can be expressed as constraints on 
the parameters of the quark matter EoS. 

Secondly, in Sec. 5, we propose a mechanism for the damp¬ 
ing of density oscillations in hybrid stars with a sharp inter¬ 
face. Such damping is important for the occurrence and satu¬ 
ration of r-modes [6,7], which if undamped can spin down a 
star in days via emission of gravitational radiation. The rele¬ 
vant observables in this case are the spin frequency, spindown 
rate, and temperature of the star. Our proposed damping mech¬ 
anism arises from the periodic conversion between quark mat¬ 
ter and nuclear matter induced by the pressure oscillations in 
the star. We find the damping grows nonlinearly with the am¬ 
plitude of the oscillation and is powerful enough to saturate an 
r-mode at very low saturation amplitude, of order 10 l0 , which 
is compatible with currently-available observations of neutron 
star spin frequencies and temperatures. 


2 The CSS parameterization 

The equation of state (EoS) of quark matter at the densities and 
temperatures relevant to neutron stars has not yet been calcu¬ 
lated from first principles. There are many models available, 
each based on some simplification of the real physics, and each 
with their own parameters. It is therefore useful to have a gen¬ 
eral parameterization of the EoS which can be used as a generic 
language for relating different models to each other and for ex¬ 
pressing experimental constraints in model-independent terms. 




2 


Mark G. Alford, Sophia Han: Characteristics of hybrid compact stars with a sharp hadron-quark interface 


We will use the “Constant-Speed-of-Sound” (CSS) parameter¬ 
ization [ 8 - 10 ], which can be viewed as the lowest-order terms 
of a Taylor expansion of the quark matter EoS about the transi¬ 
tion pressure p lrans . As illustrated in Fig. 1, one can then write 
the high-density EoS in terms of three parameters: the pressure 
/brans of the transition, the discontinuity in energy density A e at 
the transition, and the speed of sound Cq M in the high-density 
phase. For a given nuclear matter EoS £nm(p), the full CSS 
EoS is then 

£ ( P<Pt rans (1) 

^ y ^NM(/^trans) ~\~ + Cqm(P /^trans) P ^ Ptians 

where £nm(p) is the nuclear matter equation of state. 

The assumption of a density-independent speed of sound 
is accurate for a large class of models of quark matter, includ¬ 
ing many Nambu-Jona-Lasinio models [9,11-13], the pertur¬ 
bative quark matter EoS [14], and the quartic polynomial pa¬ 
rameterization [15]. For details, see Ref. [16], It is noticeable 
that models based on relativistic quarks tend to have Cq M ss 
close to 1/3, which is the value for systems with conformal 
symmetry, and it has been conjectured that there is a funda¬ 
mental bound < 7 )V] <1/3 [17], although some models violate 
that bound, e.g. [18,19] or [13] (parameterized in [9]). 



Pressure 

Fig. 1. Equation of state e(p) for dense matter. The quark matter EoS 
is specified by the transition pressure ptrans, the energy density dis¬ 
continuity Ae, and the speed of sound in quark matter Cq M (assumed 
density-independent). 


At pressures below p trans we use two examples of nuclear 
EoS: Dirac-Brueckner-Hartree-Fock (DBHF) [22] and Brueck- 
ner-Hartree-Fock (BHF) [23]. Some of the properties of BHF 
and DBHF are summarized in Table 1, where L is related to the 
derivative of the symmetry energy So with respect to density 
at the nuclear saturation density, L = 3no(dSo/dn)\„ 0 . BHF is 
a softer equation of state, with a lower value of L and lower 
pressure at a given energy density. DBHF is a stiffer EoS, with 
higher pressure at a given energy density. It yields neutron stars 
that are larger, and can reach a higher maximum mass. With 
CSS parametrization for quark matter and realistic nuclear mat¬ 
ter EoSs, we study in detail what possible forms of hybrid star 


property 

BHF, Avis 
+ UVIX TBF 

DBHF, 
Bonn A 

saturation baryon density ;io(fm _:i ) 

0.16 

0.18 

binding energy/baryon E/A (MeV) 

-15.98 

-16.15 

compressibility Kq (MeV) 

212.4 

230 

symmetry energy So (MeV) 

31.9 

34.4 

L = 3n 0 (dS 0 /dn)\ no (MeV) 

52.9 

69.4 

maximum mass of star (M 0 ) 

2.03 

2.31 

radius of the heaviest star (km) 

9.92 

11.26 

radius of M = 1.4M 0 star (km) 

11.77 

13.41 


Table 1. Calculated properties of symmetric nuclear matter for the 
BHF and DBHF nuclear equations of state used here. BHF is softer, 
DBHF is stiffer. 


mass-radius relations are and how they rely on the high-density 
phase parameter values. 

In Sec. 3 we show how the CSS parameterization is con¬ 
strained by observables such as the maximum mass M max , the 
radius of a maximum-mass star, and the radius R\ 4 of a star of 
mass 1.4M e , based on BHF and DBHF nuclear matter EoSs. 


3 Generic conditions for stable hybrid stars 

A compact star will be stable as long as the mass M of the 
star is an increasing function of the central pressure p cen t [24], 
By performing an expansion in powers of the size of the quark 
matter core, it has been shown [20,21,25] that the appearance 
of a quark matter core when p cent rises above p trans will desta¬ 
bilize the star only if the energy density discontinuity Ae is too 
large, specifically if Ae > A£ cnl where 

A £ C rit 1 3 /Jtrans 

Grans 2 2 £trans 

(This is A cr it — 1 in the notation of Ref. [20].) 

In Fig. 2, panels (b) and (c) show possible forms of M(R) 
for Ae < A£ crit , where a stable connected hybrid branch con¬ 
tinues on from the hadronic branch of the M{R) relation. Panels 
(a) and (d) show possible forms for Ae > A£ ml , where there is 
no stable connected branch, although in (d) there is a discon¬ 
nected one. In Figs. 2 (b) and (d), we illustrate how, as noted 
in Ref. [20], a second, disconnected, branch of stable hybrid 
stars can arise. The disconnected branch is a “third family” [26, 
27] of compact stars besides neutron stars and white dwarfs. In 
M{R) curves for realistic neutron star equations of state, the 
cusp that occurs where an unstable hybrid branch meets the 
hadronic branch is exaggerated in Fig. 2. In M(R) plots for re¬ 
alistic EoSs the cusp is invisibly small, covering a range in M 
of less than one part in a thousand. 


3.1 Phase diagram at fixed Cq M 

In Fig. 3 we plot a “phase diagram” for hybrid stars showing 
how the form of the M(R) relation depends on the CSS param¬ 
eters Ptrans/Grans and Ae / e lmns at fixed Cq M (see Ref. [8,16]). 
The left panel of Fig. 3 is the result of calculations for the BHF 
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Fig. 2. (Color online) Schematic form of possible mass-radius relations for hybrid stars. The thick (green) line is the hadronic branch. Thin 
solid (red) lines are stable hybrid stars; thin dashed (red) lines are unstable hybrid stars. In (a) the hybrid branch is absent. In (c) there is a 
connected branch. In (d) there is a disconnected branch. In (b) there are both types of branch. In realistic neutron starM(I?) curves the cusp that 
occurs in cases (a) and (d) is much smaller and harder to see [20,21], 


BHF + CSS(c 2 = 1) 

Htrans/Ho 




Fig. 3. (Color online) Phase diagram for hybrid star branches in the mass-radius relation of compact stars. 


nuclear matter EoS and quark matter with Cq M = 1. The right 
panel is a schematic showing the form of the mass-radius rela¬ 
tion in each region of the diagram. The regions correspond to 
different geometries of the hybrid branch displayed in Fig. 2: 
A=“Absent”, C=“Connected”, D=“Disconnected”, B=“Both” 
(connected and disconnected). 

The solid straight (red) line is Afc' m , from Eq. (2). Below the 
line in regions B and C, there is a hybrid star branch connected 
to the nuclear star branch. Above the line in regions A and D, 
there is no connected hybrid star branch. In regions B and D 
there is a disconnected hybrid star branch. 

The roughly vertical dash-dotted curve in Fig. 3 marks a 
transition where the disconnected branch of hybrid stars ap¬ 
pears/disappears. When one crosses this line from the right, go¬ 
ing from region A to D by decreasing the nuclear/quark transi¬ 
tion density, a stationary point of inflection appears in M(p cent ) 
at p cent > ptrans- If one crosses from C to B then this point of 
inflection is at higher central pressure than the existing maxi¬ 
mum in M (p C ent) ■ This produces a stationary point of inflection 


in the M(R) relation to the left of the existing maximum (if 
any). After crossing the dash-dotted line the point of inflection 
becomes a new maximum-minimum pair (the maximum be¬ 
ing further from the transition point), producing a disconnected 
branch of stable hybrid stars in regions B and D. Crossing the 
other way, by increasing the transition pressure, the maximum 
and minimum that delimit the disconnected branch merge and 
the branch disappears. 

The roughly horizontal dashed curve in Fig. 3 which sep¬ 
arates region B and C marks a transition between mass-radius 
relations with one connected hybrid star branch, and those with 
two hybrid star branches, one connected and one disconnected. 
Crossing this line from below, by increasing the energy den¬ 
sity discontinuity, a stationary point of inflection in M(p cent ) 
(or equivalently in M(R)) appears in the existing connected hy¬ 
brid branch. Crossing in to region B, this point of inflection 
becomes a new maximum-minimum pair, so the connected hy¬ 
brid branch is broken in to a smaller connected branch and a 
new disconnected branch. The maximum of the old connected 
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BHF & DBHF + CSS (c 2 = 1/3) BHF & DBHF + CSS (c 2 = 1) 




Fig. 4. (Color online) Phase diagram like Fig. 3, showing that the phase boundaries are very insensitive to changes in the nuclear matter EoS, 
but they are affected by varying the quark matter speed of sound. Note that, for a given pressure, BHF and DBHF have very different baryon 
densities. 


branch smoothly becomes the maximum of the new discon¬ 
nected branch. If one crossed the dashed line in the opposite 
direction, from B to C, the maximum closest to the transition 
point would approach the minimum and they would annihilate, 
leaving only the more distant maximum. 

Where the horizontal and vertical curves meet, the two max¬ 
ima and the minimum that are present in region B all merge to 
form a single flat maximum where the first three derivatives of 
M(R) are all zero. 


3.2 Varying Cq M and the nuclear EoS 

In Fig. 4 we show the effects of varying the speed of sound and 
the nuclear matter EoS. Fig. 4(a) is the phase diagram with 
Cqm = 1/3 (characteristic of very weakly interacting massless 
quarks), and Fig. 4 (b) is for Cq M = 1 (the maximum value con¬ 
sistent with causality). In both panels we show the phase di¬ 
agram for BHF and DBHF nuclear matter EoSs. The straight 
line is independent of Cq M and the detailed form of the nuclear 
matter EoS. The other phase boundaries, outlining the region 
where there is a disconnected hybrid branch, are remarkably 
insensitive to the details of the nuclear matter EoS but depend 
strongly on the quark matter speed of sound. For a given nu¬ 
clear matter EoS the hybrid branch structure is determined by 
/bransMrans. ^e/Grans, and Cq M , so one could make a three- 
dimensional plot with c^ M as the third axis, but this figure ade¬ 
quately illustrates the dependence on Cq M . We will now discuss 
the physics behind the shape of the phase boundaries. 


3.3 Physical understanding of the phase diagram 

The main feature of the phase diagram is that a disconnected 
branch is present when the transition density is sufficiently low. 


and the energy density discontinuity is sufficiently high. It oc¬ 
curs more readily if the speed of sound in quark matter is high. 

When a very small quark matter core is present, its greater 
density creates an additional gravitational pull on the nuclear 
mantle. If the pressure of the core can counteract the extra pull, 
the star is stable and there is a connected hybrid branch. If the 
energy density jump is too great, the extra gravitational pull is 
too strong, and the star becomes unstable when quark matter 
first appears. However, if the energy density of the core rises 
slowly enough with increasing pressure (i.e. if c 2 = dp/de is 
large enough), a larger core with a higher central pressure may 
be able to sustain the weight of the nuclear mantle above it. 
Region B, with connected and disconnected branches, is more 
complicated and we do not have an intuitive explanation for it. 

We can now understand why the vertical line marking the 
B/C and D/A boundaries moves to the right as q=, M increases. 
Since c 2 = dp/de, if </ )M is larger then the energy density of 
the core rises more slowly with increasing pressure, which min¬ 
imizes the tendency for a large core to destabilize the star via 
its gravitational attraction. Finally, we can see why that line has 
a slight negative slope: larger Ae makes the quark core heav¬ 
ier, increasing its pull on the nuclear mantle, and making the 
hybrid star more unstable against collapse. 

4 Constraining a generic high-density 
equation of state 

4.1 Maximum mass of hybrid stars 

In Fig. 5, we show how mass measurements of neutron stars 
can be expressed as constraints on the CSS parameters; the up¬ 
per plots are for DBHF and the lower plots are for BHF nu¬ 
clear matter EoS. Each panel shows dependence on /brans/Grans 
and Ae/e trans for fixed Cq M . The region in which the transi¬ 
tion to quark matter would occur below nuclear saturation den¬ 
sity («trans < «o) is excluded (hatched band at left end) because 
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DBHF (stiff) NM, Cq M = 1/3 

Btrans/Bo 


1.0 2.0 3.0 4.0 5.0 5.5 



DBHF (stiff) NM, Cq M = 1 

Btrans/flo 


1.0 2.0 3.0 4.0 5.0 5.5 



BHF (soft) NM, Cq M = 1/3 
Btrans/rio 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 



1.0 2.0 3.0 


BHF (soff) NM, Cq M = 1 

Btrans/flo 

4.0 5.0 6.0 7.0 



Fig. 5. (Color online) Contour plots showing the maximum hybrid star mass as a function of the CSS parameters of the high-density EoS. Each 
panel shows the dependence on the CSS parameters Ptrans/ e trans and 2le/£trans- The left plots are for Cq M = 1/3, and the right plots are for 
Cqm = 1- The top row is for a DHBF (stiff) nuclear matter EoS, and the bottom row is for a BHF (soft) nuclear matter EoS. The grey shaded 
region, where M m ax < 1 .95 Mq , is excluded by the observation of stars of mass M ~ 2 Mq . The hatched band at low density (where ntrans < no) 
is excluded because bulk nuclear matter would be metastable. The hatched band at high density is excluded because the transition pressure is 
above the central pressure of the heaviest stable hadronic star. 


in that region bulk nuclear matter would be metastable. There 
is also an upper limit on the transition pressure, which is the 
central pressure of the heaviest stable nuclear matter star. This 
depends on the hadronic EoS that has been assumed. 

The contours show the maximum mass of a hybrid star as a 
function of the EoS parameters. The grey shaded region, where 
7ff m ax < 1 .95 Mq, contains quark matter EoSs that are excluded 
at the one to two a level by the recent observations of stars 
of mass M ss 2M 0 [28] [29] (in Ref. [16] a slightly higher 
bound M max < 2M e was employed). For high-density EoSs 
with Cq M = 1 (right-hand plots), this region is not too large, and 
leaves a good range of transition pressures and energy density 


discontinuities that are compatible with the observation. How¬ 
ever, for high-density matter with c^ M = 1/3 (left-hand plots), 
which is the typical value in many models (see Sec. 1), the 
M max // 2M 0 constraint eliminates a large region of the CSS 
parameter space. As one would expect, the stiffer EoS gives 
rise to heavier (and larger) stars, and therefore allows a wider 
range of CSS parameters to be compatible with the 2M 0 mea¬ 
surement. The region B, where connected and disconnected hy¬ 
brid star branches can coexist, is excluded for Cq M ^ 1 /3, and 
even for larger Cq M it is only allowed if the nuclear matter EoS 
is sufficiently stiff. 
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DBHF (stiff) NM, Cq M = 1/3 
ritrans/flo 


1.0 2.0 3.0 4.0 5.0 5.5 



1.0 2.0 


DBHF (stiff) NM, c^ M = 1 

Btrans/flo 

4.0 


3.0 


5.0 5.5 



BHF (soft) NM, Cq M = 1/3 

Btrans/Ho 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 



1.0 2.0 3.0 


BHF (soft) NM, Cq M = 1 

Btrans/flo 

4.0 5.0 6.0 7.0 



Fig. 6. (Color online) Contour plots showing the radius of the maximum-mass star as a function of the CSS parameters. Dashed lines are 
for the case where this star is on the disconnected branch; for dot-dashed lines it is on the connected branch. The grey shaded region where 
M max < 1.95 Mg is excluded by the measurement of a 2M 0 star. The hatched band at low density (where ;; trans < «o) is excluded because bulk 
nuclear matter would be metastable. The hatched band at high density is excluded because the transition pressure is above the central pressure 
of the heaviest stable hadronic star. For a magnified version of the low-transition-pressure region for Cq M = 1 /3, see Fig. 8. 


In Fig. 5 the dot-dashed (red) contours are for hybrid stars 
on a connected branch, while the dashed (blue) contours are 
for disconnected branches. When crossing the near-horizontal 
boundary from region C to B the maximum mass of the con¬ 
nected branch smoothly becomes the maximum mass of the 
disconnected branch, therefore the red contour in the C region 
smoothly becomes a blue contour in the B and D regions. When 
crossing the near-vertical boundary from region C to B a new 
disconnected branch forms, so the connected branch (red dot- 
dashed) contour crosses this boundary smoothly. 

If, as predicted by many models, c^ M < 1 /3, then the exis¬ 
tence of a 2M 0 star constrains the other CSS parameters to 
two regions of parameter space: in the left-hand region the 


transition occurs at a fairly low density n tIans < 2no', in the 
right-hand allowed region the transition pressure is high, and 
the connected branch contours are, except at very low Ae, al¬ 
most vertical, corresponding to EoSs that give rise to a very 
small connected hybrid branch, and the maximum mass on this 
branch is very close to the mass of the purely-hadronic mat¬ 
ter star with p ceal = /% ans . The mass of such a purely-hadronic 
star is naturally independent of parameters that only affect the 
quark matter EoS, such as Ae and Cq M , so the contour is ver¬ 
tical. These hybrid stars have a tiny core of the high density 
phase and cover a tiny range of masses, of order 10 3 M-. or 
less (see e.g.. Fig. 5 in [8]), and so would be very rare. 
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DBHF (stiff) NM, Cq M = 1/3 
ritrans/flo 



DBHF (stiff) NM, Cq M = 1 

Btrans/flo 


1.0 2.0 3.0 4.0 5.0 5.5 



BHF (soft) NM, Cq M = 1/3 

rUrans/flo 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 



BHF (soft) NM, Cq M = 1 

Btrans/llo 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 



Fig. 7. (Color online) Contour plots similar to Fig. 6 showing the radius of a hybrid star of mass M = 1.4M 0 as a function of the CSS 
parameters. Such stars only exist in a limited region of the space of EoSs (delimited by dashed (magenta) lines). Outside that region the only 
1.4 Mq star is a hadronic star with radius 11.8 km (BHF) or 13.4 km (DBHF) (see Table 1). The grey shaded region where M max < 1.95M 0 is 
excluded by the measurement of a 2M 0 star. For a magnified version of the low-transition-pressure region for c^ M = 1 /3, see Fig. 8. 


4.2 Typical radius of hybrid stars 

4.2.1 Minimum radius 

In Fig. 6 we show contour plots of the radius of the maximum- 
mass star (on either a connected or disconnected hybrid branch) 
as a function of the CSS quark matter EoS parameters. Since 
the smallest hybrid star is typically the heaviest one, this allows 
us to infer the smallest radius that arises from a given EoS. 

The smallest stars, with radii as small as 9 km, occur when 
the high-density phase has the largest possible speed of sound 
Cqm = 1- They are disconnected branch stars arising from EoSs 
having a low transition pressure (ntrans 1$ 2no) with a fairly 
large energy density discontinuity (Ae/e tra ns > 1). The con¬ 
tours in the high-transition-pressure region are almost vertical 
because the hybrid branch is a very short extension to the nu¬ 
clear mass-radius relation, and its radius decreases with /; Irans 
in this region, because it decreases with central pressure. 


For Cq M = 1/3, the allowed low-transition-pressure region 
is disconnected from the high-transition pressure region and 
is so small that it is hard to see on this plot. By magnifying 
it (left-hand plots of Fig. 8) we see that in this region the ra¬ 
dius contours closely track the border of the allowed region 
(the M max = 1.95M 0 line) so we can say that the radius must 
be greater than 11.25 km. For a stiff hadronic EoS this mini¬ 
mum is raised to 11.4 km. These values are comparable to the 
minimum radius of about 11.8 km found in Ref. [17], which ex¬ 
plored a larger set of hadronic EoSs but did not explore the full 
CSS parameter space for the high-density EoS. If a star with ra¬ 
dius smaller than this minimum value were to be observed, we 
would have to conclude that either the transition occurs outside 
the low-density region or that Cq M is greater than 1/3. In the 
magnified figure we also show how the excluded region would 
grow if a 2.1M Q star were to be observed (long-dashed line 
for connected branch stars and short-dashed line for discon¬ 
nected branch stars). This would increase the minimum radius 
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DBHF (stiff) NM, Cq M = 1/3 
fltrans/rio 


0.5 1.0 1.5 2.0 



BHF (soft) NM, Cq M = 1/3 

Btrans/llo 


0.5 1.0 1.5 2.0 2.5 



DBHF (stiff) NM, c^ M = 1/3 
fltrans/flo 


0.5 1.0 1.5 2.0 



BHF (soft) NM, Cq M = 1/3 
fltrans/flo 


0.5 1.0 1.5 2.0 2.5 



Fig. 8 . (Color online) Magnified version of the Cq M = 1/3 plots in Figs. 6 and 7. In the two left panels, the contours are for the radius of the 
maximum mass star, which is typically the smallest star for a given EoS. In the two right panels, the contours are for R 1 . 4 , the radius of a 
1.4 Mq star. The region under and to the left of the hatched bar is probably unphysical because n trans < n o, and it was excluded (hatched band) 
in earlier figures. The grey shaded region where Af max < 1.95 Mq is excluded by the measurement of a 2 Mq star. The dashed line shows how 
that region would grow if a 2.1M 0 star were observed. 


to about 12.1 km for the soft hadronic EoS and 12.2 km for the 
stiff hadronic EoS. 


4.2.2 Radius of a 1.4M e star 

In Fig. 7 we show contours (the U-shaped lines) of typical ra¬ 
dius of a hybrid star, defined as Ri 4 , the radius of a star of mass 
1.4M 0 , as a function of the CSS parameters. The contours only 
fill the part of the CSS parameter space where there are hybrid 
stars with that mass. The dashed (magenta) lines delimit that 
region which extends only up to moderate transition pressure. 

The overall behavior is that, at fixed Ae/e tians , the typi¬ 
cal radius is large when the transition density is at its lowest. 
As the transition density rises the radius of a 1.4M 0 star de¬ 
creases at first, but then increases again. This is related to the 
previously noted fact [30] that when one fixes the speed of 


sound of quark matter and increases the bag constant (which 
increases /Arans/fimms and also varies Ae/e lvws in a correlated 
way) the resultant family of mass-radius curves all pass through 
the same small region in the M-R plane: the MIR) curves “ro¬ 
tate” counter-clockwise around this hub (see Fig. 2 of Ref. [30]). 
In our case we are varying /;> trans / £ tm ns at fixed A e / £ tra ns, so the 
hub itself also moves. At low transition density the hub is below 
1.4M 0 , so R\ 4 decreases with /Arans/Grans- At high transition 
density the hub is at a mass above 1.4M 0 so R\ 4 will increase 

with /brans / Urans ■ 

The smallest stars occur for Cq M = 1 (right-hand plots), 
where R 1.4 > 9.5km at large values of the energy density dis¬ 
continuity, and the radius rises as the discontinuity is decreased. 
This is consistent with the absolute lower bound of about 8.5 km 
[31] for the maximally compact Cq M = 1 star obeying M max > 
2M 0 . 


























Mark G. Alford, Sophia Han: Characteristics of hybrid compact stars with a sharp hadron-quark interface 


9 


DBHF + CSS (cq M = 1/3) 




DBHF + CSS (4m = 1) 


ntrans/tto~1.5 

— n t rans/n 0 =2 

— rWn 0 =2.3 

— n tran s/no=2.365 

— n t rans/n 0 = 2.7 

— n tr ans/no = 3 

■■■ n t rans/no = 3.34 


Radius of a 1.4 M 0 star (km) 


Radius of a 1.4 M Q star (km) 


DBHF + CSS (4m =1/3) 



DBHF + CSS (4m = 1) 



Radius of a 1.4 M s star (km) 


Fig. 9. (Color online) Contour plots showing how the maximum hybrid star mass and the radius of a 1.4M 0 star vary when two of the CSS 
parameters, the transition density and energy density discontinuity, are varied. Shaded regions indicate where no 1.4M 0 hybrid star exists. In 
the upper panels, EoSs below the horizontal lines (i.e. those with M max < 1.95M 0 ) are ruled out by observation, and the vertical (red) lines 
mark the radius of a 1.4M 0 purely-hadronic star. Lower panels are zoomed in on the encircled region in the upper panels, where EoSs can give 
a 1.4M 0 hybrid star on both connected branch (solid curves on the right) and disconnected branch (dashed curves on the left). 


For Cq M = 1/3 the allowed region at low transition pressure 
is small, so in the right panels of Fig. 8 we show a magnification 
of this region. We see that in the allowed (M ma x > 1.95M e and 
Htrans > «o) region there is a minimum radius 12 km for the BHF 
(soft) hadronic EoS, and about 12.25km for the DBHF (stiff) 
hadronic EoS. This minimum is attained at the lowest possible 
transition density, / i trans « hq. As the transition density rises to 
values around 2 no- the minimum radius rises to 12.4 km (BHF) 
or 13.3 km (DBHF). This is comparable to the minimum ra¬ 
dius of about 13 km found in Ref. [17], which explored a wider 
range of hadronic EoSs but assumed n t rans = 2«o- These results 
are consistent with the lower bound on R\ a for c^ M = 1/3 of 
about 11 km established in Ref. [31] (Fig. 5) using the EoS that 
yields maximally compact stars (corresponding to CSS with 
Arans = 0 and Cqm = V 3 ) obeying M max > 2M 0 . If a 1.4M 0 
star were observed to have radius below the minimum value, 
one would have to conclude that it could only arise from a CSS- 


type EoS if that EoS had Cq M >1/3. The dashed line shows 
how the excluded region would grow if a star of mass 2 . 1 M 0 
were to be observed. This would increase the minimum radius 
to about 12.7 km (BHF) or 13 km (DBHF). 


4.3 Maximum mass vs. typical radius 

Following Ref. [32], in Fig. 9 we characterize each CSS EoS by 
the maximum mass M max and the radius R\ 4 of a 1.4M 0 star 
on its mass-radius relation. We can then see which areas of the 
Ri. 4 -M max plane can be populated by hybrid stars arising from 
typical nuclear EoS combined with the CSS family of quark 
matter EoSs. Along each curve the transition density is fixed 
and we vary Ae/e trans . 

Each of the monotonic curves in Fig. 9 represents a family 
of CSS EoSs with fixed transition density, and varying A e/£ tra ns. 
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As we increase Ae/e tran s (from right to left along the curve), 
the maximum mass on a connected branch (solid curve) de¬ 
creases until it smoothly becomes the maximum mass on a 
disconnected branch (dashed curve) at the long-dashed (green) 
line, which corresponds to the nearly-horizontal phase bound¬ 
ary between region B and region C in the phase diagram Fig. 3. 
For higher transition densities the whole curve is dashed, be¬ 
cause in that case any mass-radius curve that includes a 1.4M 0 
hybrid star will always have its heaviest star on the discon¬ 
nected branch. 

The thin vertical (red) line is where the transition pressure 
has risen so high that it is equal to the central pressure of a 
nuclear matter star with a mass of 1.4M 0 (at )i tran s = 2.365 «o 
for DBHF nuclear EoS). Some of the mass contours extend to 
the right of the vertical (red) line, where a 1.4 M 0 hybrid star 
can have larger radius than the purely-hadronic star. This corre¬ 
sponds to phase transitions at low densities (« trans 2r?o) with 
a small energy density jump (Ae/e tra ns < 0.1), when at central 
pressures above transition the hybrid star M-R relation starts to 
acquire a similar behavior to that of strange stars, i.e. the ra¬ 
dius increases with the mass. At higher transition pressures, al¬ 
though on the connected branch hybrid stars are always above 
1.4M 0 , a 1.4M 0 hybrid star on the disconnected branch can 
still exist (see contours in upper panels of Fig. 9, where the 
whole curve is dashed for n trans > 2.365 «o). 

On the lower panels we show contours for EoSs that can 
give hybrid stars on both connected branch (solid curves on 
the right) and disconnected branch (dashed curves on the left), 
at transition densities very close to but below the central den¬ 
sity of a nuclear matter star with a mass of 1.4M 0 (n t rans < 
2.365 no)- From these plots we can see that with the same maxi¬ 
mum mass, radii for 1.4M 0 hybrid stars on connected/disconn¬ 
ected branches can differ by from 0.1 to 0.4km. 

Each plot in Fig. 9 contains a black horizontal line at M max = 
1.95M 0 , so all EoSs that lie below that line are observationally 
ruled out. If Cq M <1/3 (see left-hand plots), the range of radii 
for 1.4M 0 hybrid stars is limited to the vicinity of the purely- 
hadronic star with mass 1.4M 0 . This is another way of illus¬ 
trating the lower limit R > 12.25km found in Sec. 4.2.2 for the 
DBHF nuclear EoS. If c^ M = 1 (right-hand plots) the range of 
possible radii is somewhat bigger, R > 9.5km, compatible with 
the radius value tracking along the M max = 1.95M 0 boundary 
in Fig. 7, upper-right panel. 

5 Phase conversion dissipation in 
multicomponent compact stars 



t, t 2 T 

Time 


Fig. 10. (Color online) Diagram showing how the ideal boundary posi¬ 
tion (dashed line) and the real boundary position (solid blue line) vary 
in time. The ideal boundary is where the phase boundary would be 
if the phase conversion process equilibrated instantaneously, and it is 
determined by the instantaneous external pressure. The real boundary 
is always “chasing” the ideal boundary, with velocity given by Eq. (4) 
and (6) where 8z(t) is its distance from the ideal boundary. The real 
boundary coincides with the ideal boundary twice per cycle, at t = t\ 
and t = t 2 - 

sponse to the pressure oscillations in the star. This requires con¬ 
version of hadronic matter into quark matter and vice versa. If 
the finite rate of this conversion produces a phase lag between 
the pressure oscillation and the position of the interface, en¬ 
ergy will be dissipated in each cycle. We study the resultant 
damping and we find that this mechanism is powerful enough 
to saturate the r-mode at very low saturation amplitude, of or¬ 
der ICE 10 , and is therefore likely to be the dominant r-mode 
saturation mechanism in hybrid stars with a sharp interface. 

5.1 Hadron-quark conversion in a hybrid star 

In a one-dimensional, simplified model describing two phases 
in a cylinder separated by a sharp boundary moving in response 
to external oscillation in the context of Newtonian gravity, the 
expression of the corresponding energy dissipation in one cycle 
of oscillation as a function of the boundary velocity is given as 
[35] 


As noted in Sec. 1, it is very important to understand the dissi¬ 
pation mechanisms that suppress oscillations in neutron stars. 
One example is the damping of r-modes, which controls the 
spin down of neutron stars. Observational data indicates that 
the conventional damping mechanisms (bulk and shear viscos¬ 
ity) are insufficient to explain the existence of millisecond pul¬ 
sars with the observed temperatures [33,34]. 

Here we describe a dissipation mechanism that could ac¬ 
count for the observations by saturating the growth of r-modes 
at a very low amplitude [35]. The mechanism is based on the 
movement of the core-mantle interface of a hybrid star in re- 



/ f T _ dSx b (t ) 

\Jo P dr 


dr 


f x . , d5xb(r) , 

+ / A/? L sin(<»r)--dr 

Jo dr 


(3) 


assuming that the pressure in the low-density phase oscillates 
harmonically with amplitude A p b and frequency ft) around its 
equilibrium value /?, and x b is the equilibrium position of the 
boundary (the center of the star is at x = 0). Here we look at 
the case where the phases are quark matter and hadronic matter 
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and are interested in situations where phase conversion dissipa¬ 
tion becomes important in r-mode oscillations when their am¬ 
plitude is still fairly low, assuming 8p <C p in the region near 
the boundary. Therefore, we only need the EoS in a narrow 
pressure range around the critical pressure at Uansition. 

We use the calculational techniques developed by Olinto 
[36] to study the movement of the phase boundary in the sUange 
matter hypothesis scenario, but we are interested in conversion 
of quark matter to nuclear matter as well as nuclear matter to 
quark matter, since both processes occur as our burning front 
moves inwards and outwards periodically in response to an os¬ 
cillation in the pressure. The velocity of the boundary in the 
NM —> QM half cycle (f < t\ and t > ti in Fig. 10) is deter¬ 
mined by [35] 




(4) 


where = 5x;b — 5jq,| is h° w far the boundary is from its 
equilibrium position at the current pressure (see Fig. 10), and 
(q characterizes its typical length 


(»q/Xk)»Q 

2(y-l)gb£ c N rit ' 


(5) 


In Eq. (4), «n characterizes the strangeness fraction, Dq is the 
diffusion constant for flavor, Tq is the time scale of nonleptonic 
flavor-changing interactions, tjq gives the ratio of subthermal 
to suprathermal rates; in Eq. (5), iiq is the baryon number den¬ 
sity in equilibrated sUange quark matter, is the 

susceptibility with respect to K-ness evaluated at equilibrium, 
7 = iiq/h\ where «n is the baryon number density in nuclear 
matter, gb is the effective gravitational acceleration and eiL is 
the energy density at the phase boundary on nuclear matter side 
(see Sec IV in [35] for detailed discussions). 

The boundary velocity in the QM —► NM half cycle (t i < 
t < t 2 in Fig. 10) is 


dfTvb 

dr 


1 

b Q 


[ & / £ n] ' +277n<Sz/4j (6) 


where Sz = |5xb — 5x;b is how far the boundary is from its 
equilibrium position at the current pressure, with the typical 
length 


= 


(»n/Xk)»Q 
2(7-l)^b£ c N rit 


(7) 


area dS straddling the phase boundary at (R b , 0,(j)) is 

„ , „ f T , dSR b 

dW(a,8,<p) = dS(y— 1 )A/?n / cos(20 + cor)—-—dr, 

Jo dr 


( 8 ) 


where the position of the phase boundary 8R b (t) is assumed 
to move at its maximal speed [see Eqs. (4) and (6)], d5 = 




R b sin 6d6d(j) and the pressure oscillation is 

= t ?b£ c rit I SR ib I = £ c rit C b« I sin 2 0 cos 6 1 


(9) 


and all physical quantities are defined similar to those in the 
NM —► QM half cycle. 

Therefore with the periodic condition 5xb(f) = Sx b (2n/(0 + 
r), Eqs. (4) and (6) fully specify the movement of the phase 
boundary in response to the external pressure oscillation. Next 
we compute the energy dissipation in this process and see whether 
it is capable of saturating the r-mode. 


with Cb = y 75 §^^- 2R h/ R where Q and R are the rotation fre¬ 
quency and radius of the star. 

Integrating Eq. (8) over the full range of solid angle gives 
the total dissipation of the r-mode in one cycle of oscillation 
and hence the total power dissipated, /]| ls . The r-mode ampli¬ 
tude stops growing (saturates) when this equals the power in¬ 
jected via back-reaction from gravitational radiation P gI . 

As an illustrative example. Fig. 11 shows the dissipated 
power as a function of r-mode amplitude for a hybrid star ro¬ 
tating with frequency / = 600 Hz, with quark core size R b /R = 
0.56 and temperature T = 10 8 K. For the quark matter EoS we 
use the CSS parametrization introduced in Sec. 2, with n t ran s = 
4«q, Ae/e trans = 0.2, and Cq M = 1. The hadronic matter EoS is 
taken from Ref. [37], 

In the subthermal regime, the dissipated power first rises 
with the r-mode amplitude a as a 3 at very low amplitude, be¬ 
fore entering a resonant region with a maximum in P^/a 2 . 
At high amplitude in the suprathermal regime, the dissipated 
power is proportional to a 2 . The power in gravitational radi¬ 
ation from the r-mode P gr rises as a 2 at all amplitudes, and 
is also shown in Fig. 11 for this particular hybrid star. At low 
amplitude, the phase conversion dissipation is suppressed rel¬ 
ative to the gravitational radiation and therefore plays no role 
in damping the r-mode. If other damping mechanisms are too 
weak to suppress the r-mode, its amplitude will grow. How¬ 
ever, as the amplitude grows, the phase conversion dissipation 
becomes stronger, and in this example there is a saturation am¬ 
plitude a s at at which it equals the gravitational radiation, and 
the mode stops growing. 

Varying parameters such as the size of the quark matter 
core, rotation frequency, or temperature of the star will shift the 
curves in Fig. 11, and if the phase conversion dissipation is too 
weak then there will be no intersection point ( P gr will be greater 
than / J c ii s at all a) and phase conversion dissipation will not stop 
the growth of the mode. However, we can see from Fig. 11 that 
if saturation occurs, the resultant a sat is in the low-amplitude 
regime, where an analytical approach is available, and the sat¬ 
uration amplitude is extraordinarily low, of order 10 12 . This 
is typical of all model hybrid stars that we investigated. We de¬ 
rived the analytical expression for the dissipated power in the 
low-amplitude regime (dashed [black] line in Fig. 11), obtain¬ 
ing 




a 3 / 105 \ 3/2 7-1 (C 


1 £2 1 Rl 


5.2 R-mode damping and saturation amplitude 

During one cycle of an r-mode of amplitude a the energy dissi¬ 
pated in a radially oriented cylinder with an infinitesimal base 


15 \156nJ Ap N g b R 3 


where 


) 2 1 


. . _ 9 D n tin (s b £ c N rit 

m ~ b 2(7 1 t n «q(«n/Zk) 


( 10 ) 


(11) 
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Fig. 11. (Color online) (a) Dissipated power due to phase conversion P d j s (thick solid red curve) as a function of r-rnode amplitude a for a 
specific example hybrid star (see text), (b) The same quantity where the vertical axis now shows the ratio At first P d ; s is proportional 

to a 3 at very low amplitude (dashed line), then at some intermediate amplitude varies less quickly, with a maximum in Pdis/<* 2 - and finally 
changes to a 2 at higher amplitude. Also shown is gravitational radiation power P S r (thin solid blue straight line), which is proportional to a 2 at 
all amplitudes. The r-mode amplitude will stop growing when dissipation balances radiation, at the first point of intersection between the two 
curves. This defines the saturation amplitude a sat . 


For the class of models we have analyzed, in general because in 
nuclear matter diffusion is less efficient ( Dn/Dq « (9(10 2 )) 
and weak interactions take more time to proceed (T\/tq ss 
(9(10 2 )), therefore the QM -A NM transition half cycle dom¬ 
inates the dissipation. 


Eq. (10) allows us to assess how the strength of phase con¬ 
version dissipation depends on the various parameters involved. 
It is particularly sensitive to the size of the quark core, and this 
will be important when considering a whole family of hybrid 
stars with different central pressures and hence different core 
sizes. 


The results of such an investigation are shown in Fig. 12, 
where the solid (red) curve gives the numerically calculated 
saturation amplitude (a sat in Fig. 11) as a function of the size of 
the quark matter core in units of the star radius, Rb/R. To con¬ 
struct this curve we used the hadronic and quark matter EoSs 
of Fig. 11 and varied the central pressure, yielding a family 
of different star configurations. As Rb/R decreases, the dissi¬ 
pation power Pd i s decreases rapidly relative to the gravitational 
radiation F gr . The relative shift in the two corresponding curves 
in Fig. 11 leads to an upper limit on a sat when P gr is tangent to 
P dis . This corresponds to the end of the solid curve in Fig. 12 
at at the critical value of the quark core size, (Rb/R) c rit, 
below which the phase conversion mechanism cannot saturate 
the r-mode any more. 


The black dashed curve in Fig. 12 is the low-amplitude an¬ 
alytical approximation to a sat 


a, 


approx 


'2 22 tt 9 /2\ Jn (Xk ) 3 8b M2J2 r9 


3 3 -5 5 / 2 y Tn «q«n^q 


n 


rI 1 


4.2 x nr 11 y( )( TN „ 

' \ 1.5MeV ) \2x 10- 8 
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«N \ /_ 

2 no) V( 100MeV ) 2 

N x 3 


M 


\600MeV fm~ 3 J \lAM 0 


R 


(R b /R\ 
10km J \ 0.4 J 


gb 

gu 

2 

-8 


( 12 ) 


where Z>n = Ps x T- 2 , and g u is the Newtonian gravitational 
acceleration at the phase boundary when the quark core has 
uniform density e = £^ it (g u = ■ This approxima¬ 

tion is very accurate when the phase conversion damping is 
strong, but it does not capture the sudden weakening of that 
dissipation when, for example, the core radius becomes small. 

After exploring more CSS parameter space [«trans. A e/ £trans- 
Cq M ], we find that in the phase conversion mechanism, the most 
sensitive quark matter EoS parameter is the transition pressure 
(for example when the mass of the hybrid star is fixed), and the 
higher it is the larger r-mode saturation amplitude is predicted. 
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f=600Hz, T=10 8 K 



Fig. 12. (Color online) R-mode saturation amplitude (red solid curve) 
and its low-amplitude analytical approximation (black dashed curve) 
as a function of the radius of the quark matter core divided by the 
star radius R in a family of hybrid stars. For R\,/R < (Rb/^)crit ~ 0.38, 
damping is too weak to saturate the r-mode. At R^/R > 0.75 the hybrid 
star is unstable against gravitational collapse. The mass fraction of the 
core is in the range 0.12 < M core /M star < 0.68 for all configurations 
shown on the red solid curve. 

5.3 Range of validity of low-amplitude approximation 

The range of validity of Eqs. (10) and (12) is found by calculat¬ 
ing the next-to-leading (NLO) contribution, and requiring that 
it be less than a fraction e of the total dissipated power. We find 
that the approximation is valid when 


2 39 tt 5 G 2 glM 4 J 4 Q 6 ( R \ 
° 3 12 5 4 (7-l) 2 ( e N. t ) 2 V^bJ 



We see that the validity of the low-amplitude approximation is 
mainly determined by the size of the quark matter core and the 
rotation frequency of the star. 


5.4 Discussion 

We have described how phase conversion in a multicomponent 
compact star provides a mechanism for damping density oscil¬ 
lations, via the phase lag in the response of the interface be¬ 
tween components of different baryon densities to the applied 
pressure oscillation. The phase lag arises from the finite rate 
of interconversion between the phases, which limits the speed 
with which the interface can move. We studied the case where 
the two phases are separated by a sharp boundary (first-order 


phase transition) and analyzed the movement of the interface 
in the approximation of a steady state, neglecting additional 
acceleration effects and complicated hydrodynamic effects like 
turbulence. In particular, we studied the astrophy sic ally inter¬ 
esting case of the damping of r-mode oscillations [6,7] in a 
two-component star. We found that phase conversion dissipa¬ 
tion does not affect the r-mode instability region, because it 
vanishes as a 3 at low r-mode amplitude a. However, depend¬ 
ing on the values of relevant parameters, phase conversion dis¬ 
sipation can either saturate the r-mode at extremely low am¬ 
plitudes, a sat < 10 10 in the explicit example of hadron-quark 
transformation at the sharp quark-hadron interface in a hybrid 
star, or be insufficient to saturate the r-mode at all. The reason 
for this behavior stems, analogously to the bulk viscosity [38], 
from the resonant character of the dissipation, which is rela¬ 
tively strong when the time scale of the dissipation matches 
the time scale of the external oscillation (see Fig. 11). Whether 
saturation is possible depends therefore on the microscopic and 
astrophysical parameters, like in particular on the mass of the 
quark core, which should not be too small. 

Our main result is Eq. (8), which must be evaluated us¬ 
ing numerical solutions of Eqs. (4) and (6). We also give the 
low-amplitude analytic expressions for the power dissipated 
[Eq. (10)] and the saturation amplitude [Eq. (12)] which are 
valid when the dissipation is sufficiently strong, obeying Eq. (13) 
with e < 1. 

Our results have significant implications for astrophysical 
signatures of exotic high-density phases of matter, such as quark 
matter. The observed data for millisecond pulsars is not consis¬ 
tent with the minimal model of pulsars as stars made of nuclear 
matter with damping of r-modes via bulk and shear viscosity 
[34], Resolving this discrepancy requires either a new mech¬ 
anism for stabilizing r-modes, or a new mechanism for satu¬ 
rating unstable r-modes at a sat < 10 8 —10 7 [34,39,40]. Pre¬ 
viously proposed mechanisms have problems to achieve this. 
Suprathermal bulk viscosity and hydrodynamic oscillations both 
give a sat ~ 1 [41,42], The nonlinear coupling of the r-mode 
to viscously damped daughter modes could give a sat ~ 10 6 
to 10 3 [43,44]. The recently proposed vortex-fluxtube cutting 
mechanism [45] might give sufficiently small saturation ampli¬ 
tudes but is present only at sufficiently low temperatures T <C 
T c < 10 9 K, which could be exceeded by the r-mode (and/or 
accretion) heating [34]. One of the main results of this paper is 
that phase conversion dissipation can provide saturation at the 
required amplitude to explain millisecond pulsar data. 

Second, due to the extremely low r-mode saturation am¬ 
plitude of our proposed mechanism, hybrid stars would be¬ 
have very differently from neutron or strange stars. As dis¬ 
cussed in Ref. [34], if the known millisecond sources were 
hybrid stars then, for the low saturation amplitudes that we 
have found, they would have cooled out of the r-mode insta¬ 
bility region quickly (in millions of years) so that they would 
have very low temperatures by now. In contrast, in neutron 
stars r-modes would be present and would provide such strong 
heating that the temperature of observed millisecond pulsars 
would be Too ~ G(10 5 —10 6 ) K [34]. This prediction assumes a 
(so far unknown) saturation mechanism that would saturate the 
mode at a value a sat < K) x required by the pulsar data. This 
temperature is significantly higher than what standard cooling 
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estimates suggest for such old sources. The same holds for 
strange quark stars where the enhanced viscous damping can 
explain the pulsar data, but even in this case the star would spin 
down along the boundary of the corresponding stability win¬ 
dow which would keep it at similarly high temperatures. Mea¬ 
surements of or bounds on temperatures of isolated millisecond 
pulsars provide therefore a promising way to discriminate hy¬ 
brid stars. 


6 Summary and outlook 

In this paper we studied how strange quark matter could influ¬ 
ence the static properties and dynamic behaviors of hybrid stars 
in two ways. 

Firstly we explored the CSS (Constant-Speed-of-Sound) pa¬ 
rameterization of the quark matter EoS, which assumes a sharp 
transition from nuclear maftter to quark matter, and that the 
speed of sound in quark matter is independent of the pressure. 
Our study is intended to motivate the use of the CSS parame¬ 
terization as a framework in which the implications of observa¬ 
tions of neutron stars for the high-density EoS can be expressed 
and discussed in a way that is reasonably model-independent. 
We showed how mass and radius observations can be expressed 
as constraints on the three CSS parameters, and we found that 
the observation of a 2.OM 0 star already constrains the CSS pa¬ 
rameters significantly. 

If, as predicted by many models, c^ M <1/3, then there 
are two possible scenarios: a low-transition-pressure scenario, 
where the transition to the high density phase occurs at n tran s $ 
2 no- and a high-transition-pressure scenario. In the low-transiti- 
on-pressure case there are strong constraints on the radius of 
the star, as shown in Fig. 8. The radius of the maximum mass 
star (which is typically the smallest possible star) must be great¬ 
er than about 11.25 km, and the radius of a 1.4M-. star must 
be greater than about 12km [31]. The high-transition-pressure 
case corresponds to the white region on the right side of the 
left panels of Figs. 5, 6, 7, which tends to give a very small 
branch of hybrid stars with tiny quark matter cores, occurring 
in a narrow range of central pressures just above the transition 
pressure. 

If Cq M is larger than 1 /3 then a larger region of the CSS pa¬ 
rameter space becomes allowed. The right panels of Figs. 5,6,7 
show the extreme case where Cq M = 1. In this case the min¬ 
imum possible radius is around 9.0 km, and the radius of a 
1.4M e star must be greater than about 9.5 km. 

In the second part of the paper we studied a dissipation 
mechanism that is expected to occur in hybrid stars with a sharp 
boundary between the quark matter and the nuclear matter. The 
mechanism is phase conversion dissipation, and we found that 
it can saturate an unstable r-mode with an amplitude in the 
range C>( 10 12 ) to C>( 10 10 ). This is of the right order to ex¬ 
plain millisecond pulsar data. We gave an analytical form of 
the saturation amplitude, with its dependence on a set of pa¬ 
rameters such as quark core size, rotation frequency and so on. 
More measurements of temperatures would give more informa¬ 
tion about the r-mode amplitude and therefore help probe the 
compositions of hybrid stars. 

As well as the quark-hadron interface in a hybrid star, any 
first-order phase transition that leads to a sharp interface be¬ 


tween two phases with different baryon densities could lead 
to phase conversion dissipation of global pressure oscillation 
modes. For example, such dissipation could occur at the bound¬ 
ary between different phases of quark matter; for the transi¬ 
tion from nuclei to nucleon, as long as there is a difference of 
baryon densities the mechanism will operate as well. 

Our discussion was limited to the case of a sharp interface, 
which is the expected configuration if the surface tension is 
large enough. If the surface tension is small, there will instead 
be a mixed-phase region where domains of charged hadronic 
and quark matter coexist. We expect that the phase conversion 
dissipation mechanism will operate in this case too, as the do¬ 
mains expand and shrink in response to pressure oscillations. 
A similar mechanism should also be relevant for the “nuclear 
pasta” mixed phases in the inner crust of an ordinary neutron 
star. In this case in addition to the slow beta-equilibration pro¬ 
cesses there may also be slow strong interaction equilibration 
processes, whose rate is suppressed by tunneling factors for the 
transition between geometric domains of different size. This 
could further enhance the dissipation. 
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